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1  Introduction 


Background 

Much  effort  is  spent  in  improving  methods  used  to  remotely  characterize 
media  or  objects,  which  may  be  referred  to  as  media  identification.  In  general, 
media  are  identified  using  electromagnetic  or  acoustic  energy.  In  the  case  of 
electromagnetic  energy,  radar  is  used  to  perform  such  tasks.  In  particular,  a 
high-resolution  radar  is  considered  well  suited  to  perform  media  identification. 
Electromagnetic  and  acoustic  methods  are  being  used  at  the  U.S.  Army  Engi¬ 
neer  Waterways  Experiment  Station  in  many  research  programs  for  automatic 
detection  and  identification  of  subsurface  mines,  moisture  contents,  and  mater¬ 
ial  characterization. 

Target  identification  sequence  is  of  three  parts;  detection,  discrimination, 
and  recognition.  Feature  extraction  from  the  reflected  electromagnetic  energy 
is  performed  in  the  latter  for  three  reasons:  (a)  to  optimize  recognition  system 
perfomiance,  (b)  to  reduce  the  amount  of  information  to  be  processed,  and 
(c)  to  ensure  robustness  or  invariance  of  the  recogmtion  system.  For  example, 
a  set  of  features  may  actually  be  composed  of  the  fast  Fourier  transform  of  the 
reflected  energy;  target's  radar  cross  section  (RCS);  and  the  permittivity,  veloc¬ 
ity  of  propagation,  permeability,  susceptibility,  and  conductivity  profiles  of  the 
medium  in  which  the  target  exists.  While  the  target's  RCS  is  produced 
through  the  application  of  radar  methods,  application  of  inverse  scattering 
theory  is  used  to  construct  the  mentioned  profiles  (Eaves  and  Reedy  1987). 

Let  a  medium  or  target  be  a  scatterer.  An  incoming  incident  wave  impinging 
upon  the  scatterer  will  produce  a  reflected  wave  and  a  transmitted  wave.  The 
direct  scattering  problem  is  to  determine  the  reflected  and  transmitted  waves  if 
the  incident  wave  and  properties  of  the  scatterer  are  known.  In  general,  the 
inverse  scattering  problem  is  to  determine  the  properties  of  the  scatterer,  given 
the  incident,  reflected,  and  transmitted  waves.  Inverse  scattering  problems 
have  numerous  practical  applications  in  seismology,  target  identification,  sub¬ 
surface  and  ground-penetrating  radar,  geophysical  sensing,  medical  imaging, 
and  nondestructive  testing.  For  example,  in  radar  or  sonar  theory  a  known 
incident  wave  and  observed  reflected  wave  are  used  to  detect  the  properties 
and  presence  of  objects.  This  is  an  inversion  process.  Also,  in  the  operations 
of  nondestructive  testing  of  media  and  remote  sensing,  it  is  desired  to  apply 
target  or  pattern  recognition  to  identify  objects,  defects,  or  any  other  kind  of 
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targets  and  discriminate  between  them  and  clutter.  In  the  case  where  electro¬ 
magnetic  energy  is  utilized  to  serve  this  purpose,  the  electromagnetic  profiles 
mentioned  earlier  in  this  section  are  considered  important  parts  of  the  recogni¬ 
tion  process.  They  are  used  in  the  recognition  process  to  characterize  the 
different  media.  From  the  physical  point  of  view  it  should  be  noted  that  dif¬ 
ferent  bodies,  media,  or  materials  possess  different  profiles  of  electromagnetic 
properties  as  seen  from  models  of  such  properties.  Take,  for  example,  the 
simple  model  for  the  permittivity  of  a  dispersive  medium  given  in  Jackson 
(1975).  In  this  case  the  equation  of  motion  for  a  charge  is  bounded  by  a  har¬ 
monic  force  and  acted  on  by  an  electric  field  and  the  effects  of  the  magnetic 
force  are  neglected.  As  a  result  the  permittivity  is  described  as  a  function  of 
the  number  of  the  molecules  per  unit  volume,  the  number  of  the  electrons  per 
molecule,  and  a  damping  constant  Thus,  the  permittivity  profile  of  every 
material  is  unique.  Accordingly,  a  target  recognition  process  that  utilizes 
reflection,  conductivity,  peraiittivity,  and  permeability  profiles  of  media  is 
considered  a  dependable  process.  Improvement  in  the  speed  with  which  the 
results  are  generated  through  this  process  and  their  accuracy  is  logically  vital 
for  any  military  mission  such  as  the  automatic  detection  and  identification  of 
subsurface  mines.  Solutions  of  inverse  problems  are  used  for  the  generation  of 
the  electromagnetic  profiles  utilized  in  automatic  recognition  processes. 

Inverse^  problems  and  their  corresponding  solutions  have  been  formulated  in 
many  fashions.  They  are  in  the  form  of  integro-differential,  partial,  or  ordi¬ 
nary  differential  equations.  Some  of  these  problems  are  formulated  in  the  time 
domain  and  others  in  the  frequency  domain.  The  electromagnetic  profiles  are 
taken  to  be  functions  of  depth  in  the  medium,  time,  or  frequency.  For  exam¬ 
ple,  various  methods  of  reconstructing  different  profiles  through  inversion  are 
presented  in  Wolfgang-Martin,  Jordan,  and  Kay  (1981);  Kristensson  and 
Krueger  (1986);  Ladouceur  and  Jordan  (1985);  Kreider  (1989);  Bleistein  and 
Cohen  (1977);  Jordan  and  Kritikos  (1973);  Ge  (1987);  and  Cui  and  Liang 
(1993,  1994b).  Examples  of  the  application  of  iterative  numerical  methods  to 
the  solution  of  inverse  problems  are  presented  in  Devaney  (1983);  Tijhuis  and 
Van  Der  Worm  (1984);  and  Chew  and  Wang  (1990).  Consequently,  methods 
of  rating  numerical  methods  used  to  solve  a  given  inverse  problem  are  accu¬ 
racy  of  data  they  generate,  efficiency  of  time  required  by  the  method,  length  of 
computational  time,  method  implementation,  and  realization.  Improvement  of 
the  inversion  process  may  be  carried  out  in  the  formula  used.  Take  for  exam¬ 
ple  the  approximate  inversion  formula  derived  in  Ladouceur  and  Jordan 
(1985).  This  formula  implicitly  relates  the  permittivity  profile  of  an  unknown 
lossless  medium  with  the  Fourier  transform  of  the  reflection  coefficient  How¬ 
ever,  this  formula  breaks  down  when  the  discontinuity  at  the  interface  between 
the  dielectric  medium  and  free  space  is  large.  It  also  breaks  down  in  the 
deeper  regions  in  the  medium.  A  netwo±  method  is  formulated  in  Cui  and 
Liang  (1993)  to  invert  the  permittivity  profile  of  a  half-space  inhomogeneous 
medium.  This  method  is  compared  with  that  in  Ladouceur  and  Jordan  (1985) 
and  it  is  shown  that  it  can  take  into  account  possible  discontinuities  at  the  free 
space-dielectric  medium  interface.  Moreover,  as  it  is  demonstrated  in  Cui  and 
Liang  (1993),  the  netwoik  method  formulation  can  be  used  in  shallow  and 
deeper  regions  as  well;  thus  some  improvement  in  the  inversion  process  has 
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been  achieved.  This  leads  to  a  better  material  characterization  and  target  iden¬ 
tification.  Enhancement  of  inversion  processes  is  also  achieved  through 
improvements  of  numerical  methods  utilized  in  such  processes. 

Examples  of  inverse  problem  formulation  may  include  the  availability  of 
the  reflection  of  a  perpendicularly  polarized,  time-harmonic  plane  wave  with  a 
known  wave  number  incident  at  a  given  angle  on  an  inhomogeneous  region 
with  unknown  permittivity.  The  permittivity  of  the  medium  is  considered  a 
function  of  the  wave  number  used  and  depth  in  the  medium.  The  inverse 
problem  is  to  construct  the  permittivity  profile  relative  to  that  of  fiee  space. 
This  construction  is  accomplished  through  the  use  of  the  complex  reflection 
coefficient.  In  this  example  the  inverse  scattering  problem  is  formulated  in  the 
frequency  domain  and  the  medium  is  considered  dispersive  (Jordan  and 
Kritikos  1973).  In  a  second  example  a  half-space  inhomogeneous  pure  con¬ 
ducting  medium  is  considered.  A  time-harmonic  plane-polarized  electromag¬ 
netic  wave  of  a  known  wave  number  is  incident  from  the  free  space  region 
upon  the  medium  at  a  known  oblique  angle.  In  this  example  the  complex 
reflection  coefficient  is  used  to  represent  the  scattering  energy,  and  it  is  used  to 
constmct  the  conductivity  profile  of  the  medium  (Ge  1987).  A  third  example 
is  formulated  in  the  frequency  domain  where  the  considered  electromagnetic 
energy  is  scattered  from  a  one-dimensional  inhomogeneous  and  lossless  dielec¬ 
tric  slab  of  a  given  thickness  embedded  in  the  vacuum.  Using  this  scattered 
energy,  an  approximate  permittivity  profile  of  the  slab  is  generated  (Cui  and 
Liang  1994b).  Another  inverse  problem  formulation  includes  scattering  by  a 
one-dimensional  inhomogeneous  lossy  dielectric  slab  of  known  thickness 
embedded  in  the  vacuum.  The  permittivity  and  the  conductivity  of  the  slab  are 
unknown  and  assumed  to  be  real.  An  electromagnetic  pulse  of  finite  duration 
is  incident  on  the  slab,  and  the  electric  field  inside  the  slab  is  presented  in  the 
form  of  an  integral  equation  (Cui  and  Liang  1993).  Approximate  permittivity 
and  conductivity  profiles  are  generated  as  a  result  of  solving  this  inverse 
problem. 


In  practice,  situations  are  rarely  encountered  in  which  a  one-dimensional 
model  for  the  inverse  scattering  problem  can  be  used,  whereas  a  two- 
dimensional  model  proves  to  be  suitable  for  a  large  class  of  problems.  A  class 
of  solutions  targets  the  general  nonlinear  inverse  problem.  Thus  far,  no  num¬ 
erical  solutions  have  been  produced  for  these  methods  apparently  due  to  their 
inherent  instabilities.  There  have  been  hosts  of  other  solutions,  especially 
those  reported  in  connection  with  elastic  waves  in  seismic  exploration,  that  use 
optimization  and  iterative  techniques  (Santosa  et  al.  1984).  Although  these  are 
valid  methods,  their  sensitivity  to  experimental  and  numerical  noise  results  in 
inaccurate  solutions  and  poor  resolution.  Also,  a  variety  of  frequency-domain 
solutions  to  the  two-dimensional  inverse  scattering  problem  exist  based  on  the 
assumption  of  weak  scattering.  The  Bom  approximation  is  one  such  method 
which  solves  the  linear  inverse  problem.  These  approaches  are  not  valid  for 
the  more  general  cases  in  which  strong  diffraction  and  multiple  scattering 
effects  are  present  (Moghaddam  1991).  Another  class  of  solutions  targets  the 
general  nonlinear  inverse  problem  for  Schrddinger-type  equations  (Newton 
1981).  This  kind  of  approach  has  been  treated  in  connection  with  evaluation 
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of  an  optical  pulse  in  a  graded-index,  single-mode,  nonlinear  optical  fiber.  In 
this  ^proach  a  differential  equation  of  the  nonlinear  Schrddinger-type  equation 
is  obtained.  The  unknown  in  this  equation  is  the  envelope  of  a  pulse  propagat¬ 
ing  through  the  optical  fiber.  As  a  result  of  applying  the  inverse  scattering 
meth(^  an  equivalent  system  of  Gel'fand-Lavitan-Marchenko  coupled  integral 
equations  is  formulated  and  an  approximate  to  the  unknown  envelope  is  gen¬ 
erated  (Frangos  and  Frantzeskakis  1993). 

Solution  of  the  inverse  problem  is  applied  to  reduce  an  earthquake  hazard. 

It  is  used  to  better  understand  the  earthquake  processes  and  the  properties  of 
the  rock  in  the  focal  region.  In  this  case  seismographs  are  used  to  detect  and 
record  ground  motion  caused  by  the  passage  of  elastic  waves  from  earth¬ 
quakes.  The  challenge  in  this  case  is  the  determination  of  the  earthquake 
source  parameters  and  the  seismic  properties  of  the  earth  from  a  set  of  obser¬ 
vations  made  near  the  earth's  surface.  In  this  application  teleseismic  obser¬ 
vations  or  P-wave  arrival  times  from  local  earthquakes  are  inverted  for  the 
purposes  stated  above  (Santosa  et  al.  1984). 


Purpose 

This  investigation  is  a  research  in  which  efforts  are  made  to  formulate  an 
inversion  scheme.  The  scheme  should  generate  accurate  electromagnetic  pro¬ 
files.  Also,  this  scheme  should  perfonn  accurately  in  shallow  regions  and  deep 
regions  as  well.  It  should  have  a  fast  rate  of  convergence.  In  this  report  this 
scheme  is  formulated  and  used  to  generate  permittivity  and  susceptibility  pro¬ 
files.  The  scheme  is  tested  in  shallow  and  deep  regions.  Also,  effectiveness, 
accuracy,  and  rate  of  convergence  of  the  formulated  scheme  are  considered.  In 
addition,  the  scheme's  immunity  to  noise  and  its  stability  are  explored.  The 
formulation  of  the  inversion  scheme  proposed  in  this  report  is  applied  to  three 
selected  cases  to  test  its  effectiveness.  However,  the  scheme's  parametric 
study,  stability,  immunity  to  noise,  and  its  application  to  higher  dimensions 
wiU  be  conducted  in  a  subsequent  investigation.  Once  the  formulated  scheme 
is  tested  and  evaluated,  it  wiU  be  used  as  a  part  of  a  target  identification  and 
classification  process.  One  objective  in  this  process  is  to  generate  permittivity, 
velocity  of  propagation,  pemieability,  and  conductivity  profiles  effectively  and' 
accurately.  A  high-resolution  radar  is  designed  to  generate  the  electromagnetic 
energy  used  in  this  process  and  to  evaluate  the  formulated  inversion  scheme. 
This  evaluation  will  be  performed  in  the  laboratory  under  controlled  conditions 
and  followed  by  an  evaluation  in  the  field.  In  this  report,  criteria  and  require¬ 
ments  of  the  high-resolution  radar  are  addressed.  The  development  of  this 
radar  wUl  be  completed  in  subsequent  woiic. 


Scope 

Different  formulations  of  the  inverse  problem  are  given  and  discussed  in 
this  report.  The  inversion  scheme  proposed  in  this  report  is  applied  to  three 
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selected  cases.  Results  from  this  application  are  compared  with  those  gener¬ 
ated  through  the  application  of  previous  approximation  methods,  analytic^ 
deductions,  or  exact  expressions.  In  the  inversion  scheme  proposed  in  tins 
report  a  spectrd  coUocation  method  found  in  Canute  et  al.  (1988)  is  utiUzed. 
Through  this  utilization  A/*  degree  interpolating  polynomials  are  constructed  to 
approximate  functions  involved  in  the  formulation.  These  polynomials  are 
defined  by  using  Legendre-Gauss-Lobatto  points  as  the  coUocation  points  and 
Lagrange  polynomials  as  the  trial  functions.  Research  materials  related  to  the 
inverse  problem  are  given  in  Chapter  2. 

In  radar  appUcations,  high-range  resolution  can  be  achieved  by  transmitting  a 
low-peak-power,  coded  pulse  of  long  duration  and  then  compressing  it  on 
reception.  A  radar  system  that  incorporates  pulse  compression  processing 
provides  improved  detection  performarvce,  reduction  in  the  mutual  interference, 
and  an  increase  in  the  system  operational  flexibiUty  (Eaves  and  Reedy  1987). 
This  information  is  used  as  a  basis  for  the  design  and  development  of  the  lugh- 
resolution  radar  that  wiU  be  used  in  the  evaluation  of  the  formulated  inversion 
scheme.  At  the  present  time  Barker  codes  are  selected  to  form  the  coded 
transmitted  pulses.  A  brief  description  and  properties  of  these  codes  are  given 
in  this  report.  The  inverse  problem  solution  is  given  in  Chapter  2.  Research 
materials  related  to  the  high-resolution  radar  are  given  in  Chapter  3.  Conclu¬ 
sions,  recommendations,  and  comments  are  given  in  Chapter  4. 
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2  Inverse  Problem  Solution 


Formulations  of  the  Inverse  Problem 

A  general  form  of  the  solution  of  the  inverse  problem  consists  of  determin¬ 
ing  the  coefficients  of  a  partial  differential  equation  from  the  knowledge  of  the 
asymptotic  behavior  of  the  solution  (Weston  1972).  Maxwell's  equation  for 
the  electromagnetic  propagation  in  a  direction  along  the  z-axis  normal  to  a  slab 
of  varying  permittivity  e(z)^  and  conductivity  o(z)  is  in  the  foUowing  form 
(Krueger  1978): 

^zz  ~  ^0  ®  (^)  Po  ”0  (1) 


where 

E{z,t)  =  electric  intensity 
Po  =  permeability 
z  =  depth  in  the  medium 


t  =  time 

In  Equation  1  the  permeability  Pq  is  constant,  i.e.,  the  medium  is  considered 
nonmagnetic.  In  this  formulation  the  permittivity  and  conductivity  of  the  med¬ 
ium  are  functions  of  depth  only.  The  slab  is  bounded  by  free  space  where  the 
permittivity  is  Eq  and  the  conductivity  is  zero.  Using  the  transformation 


=  J  ^Po  e(-y)  ds,  H  =  x(L)  (2) 

0 


1 


For  convenience,  a  list  of  symbols  is  found  in  the  Notation  (Appendix  A). 
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where 


s  =  dummy  variable  of  integration 
//  =  the  value  of  the  variable  x  when  z  =  L 
L  =  maximum  value  of  the  variable  z 
Equation  1  is  transformed  to 

E^(x,t)  -  E^fx,i)  +  A{x)EJ^x,t)  +  B{x)Efx,f)  =  0 


where 

AW  =  -i.  — ! — ,  BW  =  --Sii 
*  ^  ® 

At  time  r  =  0  the  electric  field  impinges  upon  the  inhomogeneous  medium,  i.e. 
the  slab.  In  the  process,  a  left-moving  reflected  electromagnetic  wave  is  pro¬ 
duced.  To  model  this  physical  action,  the  electromagnetic  field  E{x,t)  is  split 
into  E  *(x,0  components  by  the  decomposition 


E\x,t)  =  1  [£(x,0  T  J EJ,x^)  ds\ 


In  Equation  4  E  '(x,t)  and  E  +(x,0  are  the  left-  and  the  right-moving  electro¬ 
magnetic  waves,  respectively.  The  relationship  between  E  (x,i)  and  E  (x,0  in 
the  inhomogeneous  medium  is  defined  through 


E  ~ix,t)  =  J  r(x,s)E  \x,t  -  s) 


In  Equation  5  r(x,f)  is  the  impulse  response  of  the  inhomogeneous  medium 
(Santosa  et  al.  1984).  The  impulse  response  satisfies  the  following  equations 

rJx,t)  -  2rXx,t)  =  -B(x)r(x,0  -  1  [A(x)  +  B{x)\r  *  Hx,t)  (6) 

**  *  2 


r(x,0)  =  -i.  [A(x)  -  B(x)] 
4 
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r{H>t)  =  0 


(8) 


where 


t 

r  *  r{x,t)  =  J  r{x,s)  r{x,t  -  s)  ds 
0 

In  application  the  impulse  response  at  =  0,  i.e.,  /-(O.O,  is  given.  Conse¬ 
quently,  the  inverse  problem  could  be  stated  as  follows:  given  the  function 
K0,0  what  are  the  permittivity  and  conductivity  profiles  of  the  inhomogeneous 
medium?  Equations  6-8  along  with  the  given  function  r{Q,t)  form  a  solution 
for  the  inverse  problem.  However,  it  should  be  noted  that  one  of  the  functions 
A{x)  or  B(x)  should  be  also  known  in  order  to  generate  the  other.  In  the  case 
when  neither  of  the  functions  AOc)  nor  B(x)  is  known,  one  is  forced  to  assume 
A(x)  or  and  solve  for  the  other.  In  application  the  function  r(0,t)  is  gener¬ 
ated  from  a  measurement  procedure.  Numerical  techniques  are  applied  to 
solved  for  the  unknown  profile.  Effects  of  parametric  changes,  noise,  rate  of 
convergence,  accuracy,  length  of  computational  time,  methods  of  implementa¬ 
tion,  and  realization  are  all  merits  on  which  the  particular  technique  is 
evaluated.  Obviously,  the  assumption  of  one  of  the  profiles  along  with  the 
complexity  of  the  process  of  selecting  one  particular  numerical  technique  is 
considered  very  important  since  the  results  from  this  stage  are  used  in  the 
target  identification  process. 

In  an  alternative  method  the  inverse  problem  is  formulated  in  terms  of  the 
frequency  domain  reflection  coefficient  of  the  medium  (Cui  and  Liang  1993). 

In  this  case  a  time-harmonic  transverse  electromagnetic  (TEM)  wave  of  wave 
number  k  is  normally  incident  on  the  interface  of  a  half-space  dielectric 
inhomogeneous  medium.  The  depth  in  the  medium  is  given  the  variable  x.  At 
the  point  x  the  frequency  domain  reflection  coefficient  of  the  medium  is  p(K,x) 
and  at  the  point  x  +  dlx  it  is  p(K,x)  +  dp(Kvr).  The  medium  is  regarded  as  a 
transmission  line  network  and  the  interface  between  the  medium  and  free  space 
is  considered  a  junction  in  this  network.  The  scattering  matrix,  S,  of  this 
transmission  line  junction  is 


where 


(9) 
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(10) 


Z  (x  +  dx)  -  Z  (x) 
Z  ix  +  dx)  +  Z  (;t) 


Z(jc)  and  Z(x+dx)  are  the  wave  impedances 
lively.  Considering  the  relationships 


Z(x)  = 


1207C 

/eW 


e(x  +  dx)  =  e(x)  + 


dE(x) 

dx 


.  dx 


at  the  points  x  and  x+dx  respec- 

(11) 

(12) 


one  can  write  Equation  10  as 


_  )/  e(x)  -  \le(x  +  dx)  ^ 
^e(x)  +  ^e(x  +  dx) 


\/e(x)  - 

y/eix)  +  1  . 

deix) 

2 

^/e(;c)  _ 

)lzix)  + 

\/e(x)  +  1 

deix) 

deix) 

4e(x) 


(13) 


Using  microwave  network  theory,  the  frequency  domain  reflection  coefficient 
of  the  medium,  pfK^r),  can  be  written  as 


p(Kvr)  = 


dT  + 


[1  -  (cfT)^]  [p(K,X)  +  £ip(K,x)]  • 
1  +  p(K,x)  dT  +  (ip(K,x)  dT 


^-j2  K'/^dx 


Neglecting  (dT)^,  (dT)'^  and  dpdT  terms.  Equation  14  reduces  to 


p(Kvr)  = 


^  p(Kj:)  -I-  dp(K,x) 
1  +  p(K,jr)  dT 


^-j2  K^z(x)dx 


(14) 


or 
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p(K^)  [1  +  p(K,x)  tffl  =  +  P(K;c) 

(dT)^  +  p(k^)  +  dp(K^)]  e 

Again  neglecting  {drf-  and  using  ^  ^\+x  yields 

p  M  +  [pM]^  dT  = 

dT  +  p(K^)  +  d^{Kjc)  -  J2k^  e(x)  dx  dT  - 

p(K^)  j2K^  e(x)  dx  -  dp(K;c)  j2K^  e(;c)  dx 

where  j  is 

Neglecting  dxdT  and  dp(K^)dx  terms  gives 

pM  +  [pM]^  dT  = 

dT  +  p(K^)  +  dp(K^)  -  p(K^)  j2K\/  e(x)  dx 

or 

rfpM  +  {1  -  [pM]^}  dT 
-  p(k^)  j2K^ej(x)  dx  =  0 


^  n  -  ip(K.  .)j2) 

dx  4e(x)  dx 

-  p(K,  X)  j2Ky/  eW  =  0 


(16) 


Equation  16  is  the  Riccati  nonlinear  differential  equation.  The  common  solu¬ 
tion  of  Equation  16  is  given  in  Ulaby,  Moore,  and  Fung  (1981  as 
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When  X  =  0,  the  ftequency  domain  reflection  coefficient  of  the  medium  is 


p(K)  = 


f  1  -  [p(K^)]^ 

J  4e(^) 


dE(s) 

ds 


-j2 


e 


(18) 


In  Equation  18  p(K)  =  p(K,0),  which  is  a  measurable  boundary  condition. 
Neglecting  the  [p(K^)]^  term  in  Equation  18  gives  the  following  equation 


Now  the  following  transformation  is  introduced; 


(19) 


then 

ds 


(20) 


(21) 


Using  this  transformation  in  Equation  19  produces  the  following  is  the  result 


p(K)  =  -J - 1_ 

0  8e(^) 


g  -jKy  ^ 

ds 


(22) 
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Equation  22  is  the  Fourier  transform  of  p(K).  It  is  an  approximate  inversion 
fonnulation  that  will  produce  accurate  results  as  long  as  the  many  assumptions 
used  in  its  derivation  are  valid.  Disregarding  the  higher  order  infinitesimals 
such  as  {dTf ,  {drf ,  dpdT  ,  dxdT,  dp(yi,K)dx,  and  the  term  [p(K^)]^  will  affect 
the  result  of  the  inversion  process  when  the  physical  meaning  of  all  or  one  of 
these  terms  is  violated.  As  a  result  a  target  identification  process  is  produced 
that  is  contaminated  with  a  high  rate  of  false  alarm. 


Pseudospectral  Legendre  approach  for  solving  Inverse  problems 

This  ^proach  for  solving  inverse  problems  is  based  on  the  idea  of  relating 
collocation  nodes  to  the  stmcture  of  orthogonal  polynomials.  Here,  polyno¬ 
mial  approximations  are  used  to  constmct  unknown  profiles.  These  construc¬ 
tions  are  made  in  terms  of  the  profiles'  values  at  Legendre-Gauss-Lobatto 
nodes. 

Let  Lf^t),  -1  <t<  1,  denote  the  Legendre  polynomial  of  order  N.  Then  the 
Legendre-Gauss-Lobatto  nodes  are  defined  by  =  -7,  =  7,  and  t  =  the 

zeros  of  Lj^ '  (t);  1  <  m  <  N  - 1  where  (t)  denotes  the  first  derivative  of 
Lf^t).  No  explicit  formula  of  the  nodes  1  <  m  <  77  -  1  is  available.  How¬ 
ever,  they  can  be  computed  numerically.  Consider  a  function  F(t)  and  let  it  be 
d^ned  over  the  closed  interval  [-1,1],  its  order  interpolating  polynomial 
F^(t)  is  defined  as  follows 


F  "(0  =  E  **  (23) 

h  =  0 


where 


N(N  -H  l)Lf^t^) 

(r^  -  1) 

0  -  h)  ' 


h  =  0,1,. 


.,N, 


(24) 


are  Lagrange  polynomials  of  order  h,  with  the  property 

=  {i  (25) 

Therefore 

The  derivative,  F(t\  of  the  function  T^r)  at  the  collocation  Mints  is  approx¬ 
imated  by  the  value  of  the  derivative,  F  ^(r),  of  the  function  /^(r)  at  the 
coUocation  points  To  find  the  first  derivative  of  F^{t)  at  the  coUocation 
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points  Equation  23  is  differentiated.  The  result  is  a  matrix  multiplication 
given  in  Gottlieb,  Hussaini,  and  Orszag  (1984)  as 


5^  bh^'hOrn) 
h  =  0 


(26) 


where 


NiN  +  1) 
4 

NiN  +  1) 

-  4 

0 


_ _  m  ^  h 

m  =  h  =  0 

m^h  =  N 
otherwise 


(27) 


and  for  the  approximation  of  the  second  derivative  of  the  function  Fit)  at  the 
collocation  points,  using  the  same  method  yields  the  following: 


2 

L^itf^ 


1 


1  <  m,h  <  N  -  1;  m  ^  h 


(28) 


^"^it^)  =  -1  -JL-,  l<h<N  -l-,m  =  h  (29) 

^  1  - 

Polynomial  interpolation  based  on  Legendre-Gauss-Lobatto  points  is  well 
behaved  compared  to  that  based  on  equally  spaced  points  (Davis  1963).  Let 
the  operation  of  interpolation  given  in  Equation  26  be  /^,  i.e.,  F^it)  =  Ij^Fif). 
Qearly  is  a  linear  operator  on  C([-l,l]),  the  Banach  space  of  continuous, 
real-valued  function  on  [-1,1],  with  the  property  =  r^.  This  space  is 
equipped  with  the  uniform  norm 

ifi-  =  _i  ITs  il'"®!’  <30) 

Since  Iff  is  a  linear  operator  with  the  property  Iff  =  If^,  then  Iff  is  a  projection 
operator  whose  range  is  the  space  of  all  polynomials  of  a  degree  less  than  or 
equal  to  N  -  1.  Furthermore,  is  a  bounded  operator  on  C([-l,l])  with 
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oo 


(31) 


>'«'  =  -I  sTs  1  E  IV')I  < 

;•  =0 

It  follows  that  (Askey  1972) 

1 

lim  i  \F^{t)  -  F{t)\P  dt^O  (32) 

>  oo  “1 


for  any  F  e  C([-l,l])  and  1  <  oo.  Using  the  LP  norm  provides 

lim  If''0) -/'(Dli,,.,  !)  =0  (33) 

N  oo 

from  which  the  following  is  deduced 

sup  |F^  I  <  CO,  V  F  e  C([-l,l])  (34) 

N  P 


Thus,  if  the  interpolation  operation  is  considered  as  a  linear  operation  from 
C([-l,l])  to  Lp(-l,l)  then,  from  the  Banach-Steinhaus  Theorem,  it  follows  that 
there  exists  a  constant  g  >  0,  which  depends  only  on  p,  such  that 

iF%<gm^  (35) 

which  is  an  upper  bound  on  the  order  interpolating  polynomial  of  the 
function  F(r). 

The  following  results  are  taken  from  Gottlieb,  Hussaini,  and  Orzag  n984) 
to  show  the  accuracy  of  approximating  Ffr)  by  F  '  ^(r)  and  F(r)  by  F  ^(t). 

a.  Theorem  1.  Let  F(r)  have  a  smooth  derivatives  for  -1  <t<  1,  and 
define 

E  f  (36) 

h  =  0  4  {dt^j 

Then  there  is  a  constant  M  independent  of  F{t)  and  N  such  that 
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^yv2p-0|^l  a  >  2p 


(37) 


b.  Theorem  2.  For  aU  FeH  s  ^  0,  there  exists  a  constant  M  inde¬ 
pendent  of  F(t)  and  N  such  that 


If  -  "l  H  S  "'vK’W-l  m„  .(.1,1).  0iq<s 


(38) 


where  a^q)  =  2q  - 1/2. 

Thus  for  a  sufficiently  smooth  function,  i.e.,  F(t),  the  rate  of  convergence  of 
F '  ^(t)  to  F  (f)  is  faster  than  any  power  of  l/N.  This  result  serves  as  a 
basis  for  the  convergence  of  the  pseudospectral  Legendre  approach  used  for 
solving  selected  inverse  problems  in  this  report.  The  next  theorem  gives 
the  error  estimate  in  the  supremum  norm. 

c.  Theorem  3.  If  to  =  -L  =  1  and  fy  ,  1<  ;  <  N  -  1  are  the  zeros  of 
Lf^t)  and  F(z)  has  no  singularities  except  a  finite  number  of  poles, 
none  of  which  lie  on  [-1,1],  and  if  for  some  n,  lF(z)/z”l-40  as 
\z\-^,  then  F%)  ->F(0  uniformly  on  [-1,1]. 

Application  of  the  Legendre-Gauss-Lobatto  collocation  method  to 
the  solution  of  inverse  problems  in  dispersive  media 

In  this  section  a  numerical  technique  is  proposed  for  solving  the  inverse 
problem  in  dispersive  scattering  theory.  The  inverse  problem  is  first  formu¬ 
lated  as  a  relationship  between  the  impulse  response  R  and  the  susceptibility 
kernel  G.  The  Legendre-Gauss-Lobatto  nodes  are  used  to  construct  the 
polynomial  interpolation  to  approximate  the  solution  of  G  for  a  given  R.  This 
method  is  efficient  and  yields  very  accurate  results.  Examples  are  included  to 
demonstrate  the  accuracy  of  this  method. 

When  a  radiation  field  interacts  with  an  object,  the  scattered  field  carries 
information  about  the  object’s  structure  and  composition.  The  inverse  problem 
is  posed  by  the  retrieval  of  material  parameters  from  measurements  of  a  scat¬ 
tered  field.  Because  of  the  complexity  of  the  scattering  interaction,  recon¬ 
structed  images  of  the  scattering  object  obtained  in  real  time  are  at  best  of 
dubious  qualitative  accuracy.  In  the  context  of  electromagnetic  (EM)  scatter¬ 
ing  when  the  illuminated  object  is  a  lossless  dielectric,  accurate  and  quantita¬ 
tive  inverse  solutions  are,  as  yet,  available  only  for  certain  types  of  profiles 
(Moghaddam  1991).  Most  solutions  are  computationally  intensive  (Jordan  and 
Ladouceur  1987)  and  are  not  attractive  for  real-time  implementation. 
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In  the  case  of  electromagnetic 
dispersive  phenomenon  lies  in  the 
1975) 


wave  propagation,  the  physical  basis  for  the 
following  constitutive  relationship  (Jackson 


D(x,t)  =  Eq 


E{x,f)  +  J  G(s)E(;c,t 
0 


-  s)  ds 


(39) 


In  Equation  39,  Eq  is  the  permittivity  of  free  space.  Equation  39  is  a  relation¬ 
ship  between  the  displacement  field  D{x,t)  at  a  point  in  a  homogeneous 
medium,  the  susceptibility  kernel  G  of  the  medium,  and  the  past  history  of  the 
electromagnetic  field  at  that  point,  E(x^)  for  -«>  <  j  <  r.  The  relationship 
between  the  susceptibility  of  the  medium  and  its  permittivity  is  given  as  (Beez- 
ley  and  Krueger  1985) 


-  1  =  J  G(r)e>'  dt  (40) 

^  0 


where  w  is  the  radial  frequency.  The  inverse  problem  considered  in  this  case 
involves  determining  the  dispersive  properties  of  a  homogeneous  medium,  i.e., 
the  susceptibility  kernel  G,  by  means  of  scattering  methods.  Here  a  direct 
computational  method  for  solving  this  inverse  problem  is  introduced.  This 
approach  is  based  on  a  spectral  collocation  method  (pseudospectral  method 
(Canuto  et  al.  1988))  in  which  degree  interpolating  polynomials  are  con¬ 
structed  to  approximate  the  susceptibility  kernel.  These  polynomials  are 
defined  by  using  Legendre-Gauss-Lobatto  points  as  the  collocation  points  and 
Lagrange  polynomials  as  the  trial  functions. 

Problem  formulation.  Consider  a  homogeneous,  isotropic,  dispersive 
medium  bounded  by  the  planes  z  =  0  and  z  =  L>  0.  An  electromagnetic  plane 
wave  is  launched  in  the  region  z  <  0.  This  impinges  normally  on  the  medium, 
giving  rise  to  a  left-moving  reflected  wave  and  a  transient  field  within  the 
medium.  Assuming  E(z,t)  denotes  a  transverse  component  of  the  electric  field 
Equation  39  gives 


£:„  +  3,2  J  G(s)E(z,t 
0 


-  s)  ds 


=  0,  0  <  z  <  L 


(41) 


0 

where  c"  =  eoPo>  Eq  Po  the  permittivity  and  permeability  of  free 
space,  respectively,  and  c  is  the  speed  of  light.  In  the  region  z  <  0,  the  field 
E(z,t)  can  be  split  into 
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E(z,t)  =  E  \z,t)  +  E  ~{z,t) 


(42) 


where  E^iz.t)  and  E  '(z,t)  are  the  incident  and  reflected  fields  respectively. 
The  relationship  between  the  incident  and  reflected  fields  is  given  in  Beezley 
and  Knieger  (1985)  as 


E  ~(z,i)  =  J  R(z,t  -  s)  E  \z,s)  ds 
0 


(43) 


where  R{z,t)  is  the  impulse  response.  Using  Equations  41-43  gives  the  follow¬ 
ing  equation  (Beezley  and  Kmeger  1985): 


2cR^  =  4/?,  +  G\t)  +  G'(0)  {2R  *R  *  R)  *  G' 
*  i2R  +  R  *  R),  0  <  z  <  L,  Q<  t<T 

(44) 

R(L,t)  =  0,  0<t<T 

(45) 

R(z,0^  =  -1  G(0\  0<z<L 

4 

(46) 

where  T  is  the  temporal  length  of  the  impulse  response,  ’ 
isk  operation  denotes  a  convolution  in  time; 

=  didt  and  the  aster- 

t 

R(z,t)  *  R(z,t)  =  J  R(z,u)  R(z,t  -  u)  du 

0 

(47) 

In  the  special  case  of  semi-infinite  medium,  0  <  z  <  o®,  /?  is  independent  of  z. 
Integrating  Equation  44  with  respect  to  t  and  using  Equation  47  produces  the 
following  equation: 

4/?i(r)  +  Git)  +  2G(0  *  Riit)  +  G(r)  *  R^it) 

*  R^it)  =  0,  0<t<T 

(48) 

In  this  case  the  solution  of  the  inverse  scattering  problem  is  to  obtain  G(r)  for 
a  given  R^  (t). 

An  approximation  of  the  susceptibility  function.  In  order  to  apply 
Legendre-Gauss-Lobatto  nodes  the  following  transformations  are  introduced: 
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^  L  (P  -H  1) 


(49) 


f  ^  r  (a  +  1) 


where  the  values  of  a  and  p  are  given  on  the  interval  [-1,1].  To  solve  Equa¬ 
tion  44  the  following  approximation  of  /?(P,a)  and  its  derivatives  are  defined: 

N  N 

/?^(P,a)=  Y:  E  /-y  <t>/ (P)  (|);  (a)  (51 

»■  =  0  j  =0 
N  N 

^p^(P>a)=  Y  E  /-y  <t>'/ (P)  f>y  (a)  (52 

I  =  0  ;■  =  0 

where  (t>i(P)  and  (t»j(a)  are  trial  functions  as  defined  in  Equation  24  and 
N  N 

(P»*^)  =  E  E  ^ij  (P)  ‘t’ /  (®) 

i=0  j=0 

The  approximation  of  G(a)  and  its  derivative  are  defined  as 


^(«)  =  E  («) 

(■  =  0 


G'>)=  E  S,ri(a) 

i  =  0 

Using  Equations  51-55,  Equation  44  reduces  to 

[2R^  +  *  R^]  +  g'^  *\l  R^  +  1  R^  *  R^ 

[T  T 

At  the  Legendre-Gauss-Lobatto  nodes  p„  ,  «  =  0,  1,  2, . N,  Equation  56 

reduces  to 
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^  E  E  '•ff  Oj  ®n)  (“)  =  y  E„  '■»/ 

L  I j  -  0  ^  J  =  0 


*  4  E  *i  (“)  2*0  E  '■»;  ♦; 

^  /  =  0  i  =  0 

+  50  E  E  "n/  J  <>.•  ('')  («  -  V)  rfV 

i  =  0  j  =  0  0 

-■is  E  *.•  s-  f  f'i  w  fy  (« -  V)  -yv 

■*  /  =  0  j  =  0  0 

N  N  N 

•=  E  E  Si  ^nk  ^nj 

■'  i  =  0  i  =  0  ;  =  0 


X  (J)',-  (a)  *  J  <t);t  (V)  <!)y  (a  -  V)  rfv 


where  v  is  a  dummy  variable.  For  a  given  R  (*l,a),  0^  i  ^  N,  Equation  57  is 
used  to  generate  a  set  of  nonlinear  equations.  This  set  of  nonlinear  equations 

is  solved;  as  a  result  the  coefficients  g; .  /  =  0,  1 . .  N  are  generated,  and 

using  Equation  54,  the  solution  of  the  inverse  problem  can  be  calculated.  In 
the  special  case  of  semi-infinite  medium,  let  Q(t)  =  /?i(r)  *  Riif),  and  q  =  ¥iT 
(0  +  1)  where  the  value  of  9  is  given  on  the  interval  [-1,1].  Using  Equation 
48,  the  following  is  the  result 


fT  \  IT 

4/?i  i.  (a  +  1)  +  G  ±  (a  +  1)  +T 
2  1 2 


(2l/7)-l  \  r  T 

J  G  1(9  *  1)  Ri  t  - -0  *  1) 

-1  J  K 

I  r  G  I  (0  + 1)  Q  t  -  4-(e  1) 


J  r  -  ^  (0  +  1)  dO 


de  =  0 


At  a  =  4  ,  A:  =  0,  1,  2, . .  N, 
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AR 


N 


■TT  (^k  *  *  Sk*"^ 


8 


m 


to  N{N*\)Lf^{t0 


A 

x/ 


-1  L 


(02  -  1)  L'„  (0) 


0  -  t 


m 


R 


*2^ 


(>i  -  6) 


^  «*  -  9) 


d0  =  0 


(59) 


Equation  59  is  an  algebraic  equation  that  can  be  solved  for  the  unknowns  g- ,  i 
=  0,  1.  2,..,  N. 


Numerical  examples.  This  Legendre-Gauss-Lobatto  nodes  method  is 
applied  to  solve  two  examples.  In  both  examples  the  impulse  response  is 
known,  whereas  G(t)  needs  to  be  estimated.  For  the  purpose  of  this  demon¬ 
stration  the  direct  problem  had  to  be  solved  in  order  to  generate  R(0,t).  This 
function  then  is  used  to  solve  the  inverse  problem  and  generate  the  correspond¬ 
ing  G(t). 


a.  Example  1.  A  two-resonance  model  for  the  electron  contribution 
to  the  permittivity  is  used  in  this  example.  Thus  G  is  given  by 
(Beezley  and  Krueger  1985) 


G(r)  =  e  sin  (1.6w)  +  0.5e  "^•5'  sin  (6n)  (60) 

The  depth  of  the  medium  is  taken  as  0.8  and  the  chosen  time  T=6. 
The  approximate  value  of  the  unknown  function,  i.e.,  G  is  calcu¬ 
lated  for  N  =  4.  In  Table  1,  the  values  of  the  approximate  and  exact 
susceptibility  fonction  at  the  Legendre-Gauss-Lobatto  nodes  are 
reported. 


Table  1 

Results  of  Example  1 

(n 

-1.000000 

0.0014 

0 

-0.654654 

-0.4799 

-0.4808 

0 

0.3228 

0.3226 

0.654654 

-0.0992 

-0.0996 

1.000000 

-0.2855 

-0.2864 
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b.  Example  2.  In  this  example  the  medium  is  semi-mfinite;  the  impulse 
response  is  taken  as  =  t,  and  the  time  is  taken  as  T  =  10.  The 
exact  solution  for  G(0  is  given  as 

(< 

Gif)  =  -2  (sin  t  +  t  cos  t) 


The  approximate  and  exact  susceptibility  functions  at  the  Legendre- 
Gauss-Lobatto  nodes  for  N  =  8  are  given  in  Table  2. 


Table  2 

Results  of  Example  2 

OinSS 

G(»n) 

-1.000000 

0.0012 

0 

-0.899757 

-1.8410 

-1.8401 

-0.677186 

-1.8495 

-1.8584 

-0.363117 

6.4341 

6.4486 

0 

-0.9205 

-0.9187 

0.363117 

-12.7680 

-12.7597 

0.677186 

6.7687 

6.7832 

0.899757 

19.1291 

19.0935 

1.000000 

17.8995 

17.8695 

Application  of  the  Legendre-Gauss-Lobatto  collocation  method  to 
Gei'fand-Levitan  integral  equation 

One  way  to  solve  a  one-dimensional  inverse  scattering  problem  is  through 
the  Gei'fand-Levitan  integral  equation.  The  closed-form  solution  of  this  inte¬ 
gral  equation  is  not  easy  to  obtain.  In  this  section  a  numerical  method  for 
solving  this  integral  equation  is  proposed,  hence  an  alternative  method  of  con¬ 
struction  of  arbitrary  dielectric  profiles.  Two  examples  are  given  to  demon¬ 
strate  the  accuracy  of  this  method.  Pioneering  work  in  this  subject  was  per¬ 
formed  by  Gel'fand  and  Levitan  (1955),  Marchenko  (1955),  and  Kay  (1955, 
1960).  These  authors  determined  the  scattering  potential  associated  with  the 
time-independent  Schrddinger  equation  from  scattering  data,  such  as  reflection 
and  transmission  coefficients.  The  Gei'fand-Levitan  integral  equation  is  a 
formulation  of  an  inverse  problem.  However,  explicit  solutions  of  this  equa¬ 
tion  are  seldom  easy  to  obtain  (Kritikos,  laggard,  and  Ge  1982).  Solutions  for 
the  Gei'fand-Levitan  integral  equation  have  been  constructed  by  Kritikos, 
laggard,  and  Ge  (1982),  Jordan  and  Ladouceur  (1987),  Kay  (1960),  and  Reilly 
and  Jordan  (1981).  Recent  approximate  methods  of  solving  inverse  problems 
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arc  given  in  Cui  and  Liang  (1994a,  1995),  Razzaghi  and  Ahmad  (1995),  and 
Kincanon  (1993). 

Problem  formulation.  The  Gel'fand-Levitan  approach  is  based  on  the 
solution  of  the  following  integral  equation  (Kritikos,  laggard,  and  Ge  1982): 


R(z  +  >-)  +  K{z,y)  +  J  A"  (z,z')  R  (y  *  z')  dz'  =  0 

^  (62) 
z  >  0 

-z  <  y  <  z 

where  K  (z,y)  is  the  kernel  and  A  (z  +  >>)  is  the  transformed  reflection  coeffi¬ 
cient  The  dielectric  permittivity,  e(z),  is  related  to  the  potential  function,  v(z), 
through  the  following  relationship  (Kritikos,  laggard,  and  Ge  (1982)): 

e(z)  =  1  -  v(z)/k^  (63) 

The  solution  for  the  potential  function  produces  the  following  equalities 
(Gel’fand  and  Levitan  1955;  Kay  I960): 


vfzl  -  p-dK(,z^)/dz,  z  >  0 
“  \0,  z  <  0 

(64) 

[9^  K{z,y)/dz^]  -  [02  Kiz,y)/dy^] 

(65) 

-  [2K{z,y)  dKiz,z)/dz]  =  0 

R(z)  =  0,  z  <  0 

(66) 

The  aim  is  to  solve  for  the  permittivity  e(z).  Substituting  the  solution  of  Equa¬ 
tion  65  in  Equation  64  and  making  use  of  Equation  63  an  expression  for  e(z) 
is  produced. 

An  approximation  of  the  dielectric  profile.  In  order  to  apply  Legendre- 
Gauss-Lobatto  nodes  the  following  transformations  are  introduced 

z  =  .1  (a  +  1)  -  ^  (a-1)  (67) 
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where  zy,  Z2,  yj  and  y2  are  the  limits  of  z  and  y.  The  values  of  a  and  p  are  in 
the  interval  [-1,1].  Define  the  approximation  of  K(z,y)  and  its  derivatives  as 
follows: 


N  N 

K^ia,^)=  5^  X)  fly '!),<«)  (|)y  (P) 

1  =  0  j  =  0 


(69) 


Ka  (««P)  =  E  E  ^ij 

i  =  0  j  =  0 


(70) 


N  N 

<  (“.P)  =  E  E  Oij  ■I’/  <“)  ®> 

i  =  0  j  =  0 

Using  Equations  67  and  68,  Equation  65  takes  the  following  form: 


(71) 


^2^2  a2/i:(a,p)  _  (ly-  d^Kja,^) 

df 


yh 


da^ 


yh 


(72) 


-  2/i:(a,p) 


dK  (a,a)  _  Q 


9a 


where  L  =  Z2  -  zj  =  y2  ~  yi-  Using  Equations  69-71,  Equation  72  is  reduced  to 
N  N  N  N 

E  E  “.J  f /"  (“)  ‘tj  ®  -  E  E 

;  =  0  ;•  =  0  i  =  0  j  =  0 


N  N 

a.  <t);  (a)  if/'  (^)  -LK''  (a,P)  £  ^ 


1  =  0  j  =  0 


fly  [ifi  (a)  ifj  (a)]'  =  0 


(73) 


At  the  Legendre-Gauss-Lx)batto  nodes  a  =  and  p  =  P^j,  Equation  73 
reduces  to 
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L  a, 


N  N 

E  (a«)  -  E  (PJ 

1  =  0  y  =  0 

N 

E  S- '{’/  E  ««  <(>/'  (a„) 


nm 


U  =  o 


i  =  0 


(74) 


The  unknowns  in  Equation  74  are  the  coefficients  a-;  /,  y  =  0,  1,  2,...,  M  In 
order  to  solve  for  the  unknown  coefficients,  the  following  conditions  will  be 
used.  Differentiating  Equation  62  gives 


dR(z  *y)  ^  dK(z,y)  , 

—h  ~w~ 

j  K(z,z')  dz'  +  K{z-y)  /?(0)  =  0 

iy  ^ 

and 


(75) 


dy^  dy^ 

^  R(y  +  z^)  ,  /  j, ,  . 

±dz'  +K  (z-y) 


dy‘ 


dR  (y  +  z^) 


dy 


Z  =  -y 


+  (^~y)  R  (())  =  0 

dy 


(76) 


Let  z  =  -y,  Equations  75  and  76  reduce  to 


dR(z  +  y) 


dy 


^  dKiz,y) 


z  =  -y 


dy 


+  K(z,z)  R  (0)  =  0 


z=  -y 


(77) 


and 
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Equations  77  and  78  are  two  conditions  under  which  Equation  74  is  solved. 
As  a  special  case  when  R(0)  =  0 


Also  from  Equation  62 
K  (z-z)  =  0 


(81) 


As  n,m  =  0,  1,  2,  3,...,  N,  Equation  74  reduces  to  a  set  of  nonlinear  algebraic 
equations.  Solving  this  set  of  equations  results  in  values  for  the  unknown 
coefficients.  The  initial  value  of  the  potential  function  is  deduced  from  Equa¬ 
tions  62  and  64.  Its  final  expression  is 


Numerical  examples.  Two  examples  of  inverse  problems  are  selected  and 
solved  to  test  the  numerical  scheme  proposed  in  this  section.  The  first  exam¬ 
ple  was  considered  in  Kritikos,  laggard,  and  Ge  (1982).  It  was  solved  numeri¬ 
cally  through  an  iterative  method.  The  second  example  was  considered  in 
Moses  and  de  Ridder  (1963)  and  solved  analytically.  Given  the  complex 
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reflection  coefficient,  r(K),  the  requirement  in  both  examples  is  to  reconstruct 
the  dielectric  profile  as  a  function  of  depth,  i.e.,  e('z).  As  it  is  seen  from 
Equation  63,  this  requirement  is  translated  to  reconstructing  the  potential  func¬ 
tion  vf'z). 


a.  Example  1 .  As  a  result  of  a  test,  the  reflection  coefficient  is  given  as 


r(K) 


-^1  k2 

(K  -  (k  - 


(83) 


where 


^  (1  -J) 


and 


h 


-(1  *j) 


The  transform  of  /<k)  is  (Kritikos,  Jaggard,  and  Ge  1982) 


(1  *  })  -z  (1  -  f) 


Riz)  = 


-  -J 


yll 


(84) 


This  example  is  treated  in  the  area  bounded  by  0  <  z  <  0.5  and  0  <  y  <  0.5. 
The  exact  and  approximate  values  of  the  function  vfz)  are  given  in  Table  3. 


Table  3 

Values  of  the  Function  v(z).  Example  1 

Z 

Approximate  v(z) 

Exact  v(z) 

0.00000 

4.00000 

4.00000 

0.08634 

3.17540 

3.17686 

0.25000 

2.18449 

2.18328 

0.41366 

1.59203 

1.59220 

0.50000 

1.37141 

1.37258 

b.  Example  2.  The  reflection  coefficient  from  testing  a  certain  dielectric 
medium  is  given  as 
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(85) 


r  (K)  = 


-^1  k 

(K  -  k{)  (K  -  k2) 


where 

ki=J 
and 
h  = 

The  transform  of  Kk)  is  (Moses  and  deRidder  1963) 
r(z)  =  Sink  (z) 


This  example  is  also  treated  in  the  area  bounded  by  0  <  z  <  0.5  and 
0  <y  ^  0.5.  The  exact  and  approximate  values  of  the  function  v(z)  are  given 

in  Table  4. 


Table  4 

Values  of  the  Function  v(z) 

z 

Approximate  v(z) 

Exact  v(z) 

0.00000 

-4.00000 

-4.00000 

0.08634 

-3.94063 

-3.94096 

0.25000 

-3.53830 

-3.53890 

0.41366 

-2.89300 

-2.89205 

0.50000 

-2.51071 

-2.51716 

Application  of  the  Legendre-Gauss-Lobatto  collocation  method  to 
Zakharov-Shabat  equations 

This  is  a  case  in  which  a  numerical  method  is  proposed  to  solve  the  inverse 
scattering  problem  associated  with  two  Gel'fand-Levitan-Marchenko  integral 
equations.  This  inverse  problem  is  represented  by  two  coupled  parti^ 
differential  equations.  The  solution  of  these  partial  differential  equations  is 
two  kernel  functions.  These  kernels  are  also  related  to  each  other  through  the 
Gel'fand-Levitan-Marchenko  integral  equations.  Legendre-Gauss-Lobatto 
nodes  are  used  to  construct  the  polynomial  interpolation  to  approximate  the 
solution  of  the  kernel  functions  related  to  the  scattering  potential.  Zakharov- 
Shabat  equations  are  used  to  formulate  inverse  scattering  problems  in  Frangos 
and  Jaggard  (1991,  1992).  Gel'fand-Levitan-Marchenko  coupled  integral 
equations  are  associated  with  the  2^akharov-Shabat  equations  (Lamb  1980). 

Applications  of  the  solution  of  these  types  of  inverse  scattering  problems  are 
presented  in  Song  and  Shin  (1985). 
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Again,  this  numerical  approach  is  based  on  a  spectral  collocation  method 
given  in  Canuto  et  al.  (1988).  In  this  numerical  approach,  A/*  degree  inter¬ 
polating  polynomials  are  constructed  to  approximate  the  kernels,  which  in  turn 
are  related  to  the  scattering  potential.  These  polynomials  are  defined  through 
the  use  of  Legendie-Gauss-Lobatto  points  as  collocation  points  and  Lagrange 
polynomials  as  trial  functions. 

Problem  formulation.  The  linear  two-component  Zakharov-Shabat  coup¬ 
led  equations  are  given  in  Lamb  (1980)  as 

<A’i(:r,K) 

— - -  +  /KVi(x,K)  =  q{x)\2{x,K)  (87) 


— -  ncv2(^-’^)  =  -‘7  V)Vi(x,K) 


(88) 


where 

Vi(x,k),V2(x,k)  =  modes  of  propagation 
qix)  =  potential  function 

and  the  ^terisk  denotes  conjugation.  In  this  case  x  represents  depth  in  the 
medium  in  which  the  propagation  is  taking  place.  As  a  solution  of  this  inverse 
problem,  given  the  reflection  coefficient,  r(K),it  is  desired  to  find  the  function 
q{x).  r(K)  is  the  ratio  of  the  reflected  and  incident  electric  fields.  It  is  shown 
in  Lamb  (1980)  that  in  the  region  x  >  y  Equations  87  and  88  are  related  to  two 
Gel'fand-Levitan-Marchenko  coupled  integral  equations.  These  integral  equa¬ 
tions  are 


X 

>llV.y)  +  J  A2{x,t)R(y  +  t)  dt=0  (89) 

-y 


X 

-  A2(x,y)  +  R(x  +y)  +  j  Ai(x,t)R(y  +  t)  dt  =  0  (90) 

-y 


where  Ai(y)  and  A2(.x,y)  are  the  kernels  and  Rix)  is  the  Fourier  transform  of 
the  reflection  coefficient  r(K).  In  region  x>y  the  relationships  between  the  two 
kernels  >li(x,y)  and  A2{x,y)  are  given  in  Lamb  (1980)  as 
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— i -  +  — i - =  q(x)A2(x,y) 

dy  dx 

dA2(x,y)  dA2(.x,y) 

— —  =  q  (A:)Ai(;c,y) 

where 

q\x)  =  -2A2(xa) 


and 


\q(,xT  =  -2 


dAi  (a:^) 
dx 


From  Equations  89  and  90,  in  the  limit  y-^  -x 


A^  (x,-x)  =  0 


A2(x-x)  =  R(0) 

Differentiating  Equations  89  and  90  and  letting  y-*  ‘X 
dA^(;)c,y) 


+  y42(jC,JC)/?(0)  =  0 


(91) 

(92) 


(93) 

(94) 


(95) 


dA2ix,y) 

By 


y  = 


+ 


-X 


dR(x  +  y) 
— 


y  =  -X 


+  Aj(jc,x)/?(0)  =  0 


(96) 


Equations  93-96  are  conditions  under  which  Equations  91  and  92  are  solved 
and  solutions  for  the  two  kernels  Aj(a:,>’)  and  A2ix,y)  are  generated.  Conse¬ 
quently,  a  solution  for  the  potential  function  q(x)  is  produced. 
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An  approximation  of  the  potential  function.  In  order  to  be  able  to  apply 
Legendre-Gauss-Lobatto  nodes  and  solve  the  inverse  problem,  consider  the 
following  transformations: 


oc  =  _  (1  +  a) 


>’  =  y  (1  ^  P) 


where  L  is  the  range  of  x  The  values  of  a  and  p  are  in  the  interval  [-1,1]. 
Define  the  approximations  of  the  kernels  Aj(x,y)  and  A2(,x,y)  as 


N  N 


Ai  (a,p)=  Z  Z  Yy<t),  (a)(t)-(P) 
,•  =  0  y  =  0  ^ 


N  N 


A2  (a,p)  =  Z  Z  ([,.  (a)  ([)-  (p)  (lOO: 

,-  =  0  y  =  0  ■' 

where  y  and  C  are  two  coefficients  and  the  function  ^  is  given  in  Equation  24. 
From  the  definition  of  q  (x)  (following  Equation  92) 


N  N 


^V)=-2  Z  Z  Ty  <t»i  (a)  <t);  (a) 
/  =  0  y  =  0  •' 


The  differentiations  of  Equation  99  with  respect  to  a  and  p  are 


N  N 


Ai(a,P)„=  Z  Z  Yy  <!>/' (a)  (j);  (P) 

/  =  0  y  =  0 


N  N 


Aj(a,P)p=  Z  Z  Yy  <t>/ (a)  (|)/ (p) 

/•  =  0  y  =  0  ■' 

Equations  similar  to  Equations  102  and  103  can  be  generated  from 
Equation  100. 


Application  of  Equations  99  and  100  and  their  corresponding  derivatives  to 
Equations  91  and  92  results  in  the  following  equations 
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N  N 


N  N 


E  E  Y//  <(>,•  (a)  <(>/  (P)  +  2  E  Yy  't>,'  («)  <t»/P) 
,=oy=o  ^  ^  i=oj=o  ' 

N  N  ^ 

=  -L  E  E  C;/  <(>/  (“)  4»/  (“) 

,  =  oy  =  o  •' 


.  z  2  Cy  «t>.- (a)  «t>y  (p) 

«  =  0  ;•  =  0 


N  N 


N  W 


E  E  Cy  <t»,-  (a)  <1>'/  (P)  -  2:  E  Cy  «)>/  (a)  ifj  (P) 
.•  =  oy  =  o  ^  /  =  0;  =  0  ' 

=  -L  E  E  Cy  <1>/  (a)  (“) 

;  =  0  y  =  0 


N  N 

2  2  Yy  <!>/  (a)  (t),-  (P) 

(•  =  0  y  =  0 


Considering  the  case  where  R(x)  is  real,  at  the  Legendre-Gauss-Lobatto  nodes 
a  =  a„  and  P  =  P„  Equations  104  and  105  reduce  to 

T  I  (106) 

y  =  0  1  =  0 

2  c„j  (P;„)  -  2  c,„,  (()/  (a„)  =  -  L  c„„  .  (107) 

y=o  «=o 

where  a  and  c  are  the  real  parts  of  the  coefficients  y  and  C-  As  /i,  m  =  0,  1,  2, 
N  two  sets  of  algebraic  equations  are  generated.  Solution  of  these  alge¬ 
braic  equations  results  in  values  for  the  coefficients  a  and  c.  Consequently, 
from  Equation  101  the  potential  function  q(x)  is  known  at  the  Legendre-Gauss- 
Lobatto  nodes. 

A  numerical  example.  Let  the  reflection  coefficient  be  given  as 

riK)  =  i  -  J_  (10® 

2  (K  -J) 


The  Fourier  transform  of  /"(k)  is 
Rix)  =  le-^ 
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The  exact  expression  for  the  corresponding  potential  function  is  given  in  Fran- 
gos  and  Jaggard  (1991)  as 


q{x)  = 


2s  (1  -  s) 

{s  -  1)2  6-2^  +  I  e'^ 
4 


(110) 


where  s  =  1.11803399.  The  Legendre-Gauss-Lobatto  method  was  used  to 
solve  this  inverse  scattering  problem  and  generate  the  corresponding  potential 
function  numerically.  The  computations  were  carried  out  for  AT  =  4  and  =  6, 
Results  from  these  computations  are  presented  in  Tables  5  and  6. 


The  first  column  in  Tables  5  and  6  shows  the  Legendre-Gauss-Lobatto  nodes. 
The  second  column  shows  the  values  of  the  computed  potential  function,  and 
the  third  column  shows  the  corresponding  exact  values. 


Table  5 

Estimated  and  Exact  Values  of  q(a)  with  Af  =  4 

“n 

0  ^  n^4 

-1.000000 

-1.0000000 

-1.0000000 

-0.654654 

-0.4776921 

-0.4820000 

0.000000 

-0.1124612 

-0.1127600 

0.654654 

-0.0257691 

-0.0261026 

1.000000 

-0.0119132 

-0.0120594 

Table  6 

Estimated  and  Exact  Values  of  q(a)  with  N=6 

“#i 

0 

-1.00000000 

-1.00000000 

-1.00000000 

-0.83022390 

-0.70289322 

-0.70388000 

-0.46884879 

-0.32025671 

-0.32025800 

0.00000000 

-0.11274966 

-0.11276200 

0.46884879 

-0.03953648 

-0.03954590 

0.83022390 

-0.01758667 

-0.01762760 

1.00000000 

-0.01207034 

-0.01205940 
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Effect  of  Noise  on  the  Pseudospectral  Legendre 
Method  Applied  for  Solving  Inverse  Problems 

A  criterion  for  the  effect  of  noise  on  the  Legendre-Gauss-Lobatto  pseudo¬ 
spectral  method  is  formulated  in  this  section.  For  the  purpose  of  demonstra¬ 
tion  a  particular  inverse  problem  is  utilized.  Consider  the  inverse  problem  of  a 
semi-infinite  scattering  medium  in  the  region  z  >  0,  in  which  the  velocity 
profile  c(z)  is  arbitrary,  with  c(z)  constant  in  z  <  0.  The  gradient-type  interface 
at  z  =  0  is  characterized  by  a  jump  in  the  first  derivative  of  c(z),  i.e.,  c^.  In 
addition,  c(z)  itself  is  continuous  everywhere,  while  is  continuous  every¬ 
where  except  at  the  interface.  The  form  of  the  wave  equation  used  here  is 
(Kieider  1989) 

M2z(z,0  -  c~Hz)Ut,(z,f)  =  0 


Introducing  the  travel  time  coordinate  x 
z 

x(z)  =  J  dale  (a) 

0 


Equation  111  takes  the  following  form 

v^(x,0  -  Vflx,t)  -  b(x)v^  (x,t)  =  0 


(112) 


(113) 


where 


v(x,0  =  u(z,t) 


and 


b(x)  =  cj^z) 


Given  some  knowledge  of  the  electromagnetic  field,  the  solution  of  this 
inverse  problem  is  to  reconstruct  the  velocity  profile,  i.e.,  c(z),  starting  from 
the  interface  up  to  a  specified  distance  into  the  medium,  i.e.,  Z.  The  strategy 
is  to  let  a  known  right-moving  incident  wave  propagate  from  z  =  -«»  so  that  it 
impinges  upon  the  inhomogeneous  medium  beginning  at  time  r  =  0.  The 
impulse  response  of  the  medium,  r  (x,t),  satisfies  the  following  equation 
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(114) 


rp:,t)  -  2r,  (x,t)  =  1  b  (x)  r  *  r 


with  initial  condition 

r  (j:,0)  =  ^  6  (x) 
4 


(115) 


where  the  asterisk  denotes  a  convolution  in  time.  In  addition  the  following 
integral  is  defined 


Z 

L  =  jda/c  (a)  (116) 

0 


For  a  given  r(0,t),  one  can  apply  the  formulated  Legendre-Gauss-Lobatto 
collocation  method  and  solve  this  inverse  problem.  The  question  is  how  this 
method  will  behave  when  the  function  r(0,t)  is  contaminated  with  noise.  To 
answer  this  question,  consider  the  following  criteria; 

a.  Let  the  noise  be  additive.  Accordingly  the  following  equation 

A0,t)  =  r(0,t)  +  n(t)  (117) 


represents  the  contaminated  impulse  response  where  n(t)  is  the  additive 

noise. 

b.  n(t)  may  take  many  statistical  shapes.  Distribution  functions  may  be 
applied  such  as  the  Uniform,  Arcsin,  Beta,  Cauchy,  Chi-Square,  Erlang, 
Exponential,  Gamma,  Gaussian,  Laplace,  Log-Normal,  Rayleigh,  Rice, ' 
or  WeibuU  to  represent  the  noise. 

c.  In  every  case,  the  contaminated  impulse  response  given  in  Equation  117 
is  used  as  an  input  to  the  formulated  Legendre-Gauss-Lobatto  col¬ 
location  method,  and  the  corresponding  velocity  profile  is  generated. 

d.  Data  generated  from  step  c  can  be  used  to  make  comparisons  and  study 
the  behavior  of  the  formulated  Legendre-Gauss-Lobatto  collocation 
method  in  the  presence  of  the  different  types  of  noise.  The  presence  of 
one  type  should  be  exclusive. 

e.  Represent  the  root  mean  square  of  the  signal-to-noise  ratio  as 
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where  SD  is  the  standard  of  deviation  and 
1 

r  =  —  f  r  ds 

2  J 


(118) 


(119) 
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3  Radar  Employment 


Some  of  the  investigations  described  in  this  Chapter  are  laboratory  pro¬ 
cedures  that  will  be  discussed  in  a  subsequent  report.  However,  such  proce¬ 
dures  are  mentioned  here.  Thus  far  in  this  report  a  formulated  scheme  for 
solving  inverse  problems  is  proposed.  In  this  Chapter  the  design  and  criteria 
of  a  high-resolution  radar  are  discussed.  This  radar  is  designed  with  the  con¬ 
cept  of  its  employment  for  the  evaluation  of  inversion  schemes.  As  discussed 
in  the  preceding  Chapter,  the  input  to  the  inversion  processes  that  is  selected 
for  the  application  of  the  formulated  Legendre-Gauss-Lobatto  collocation 
method  is  the  impulse  response  or  complex  reflection  coefficient  of  the  tested 
medium. 

A  network  analyzer  that  is  configured  as  a  radar  system  will  be  used  to 
generate  the  impulse  response  of  a  medium  under  testing.  Characterization  of 
the  system’s  configurations  is  not  within  the  scope  of  this  report 

A  high-resolution  radar  is  designed  to  generate  the  complex  reflection  coef¬ 
ficient  in  the  laboratory.  This  radar  wiU  be  developed  before  the  evaluation 
and  testing  of  different  formulated  inversion  schemes.  This  radar  generates  a 
specific  carrier  wave  form.  Its  transmitted  signal  is  composed  of  a  coded 
version  of  this  wave  form.  Barker  codes  are  used  in  the  encoding  process. 

For  convenience,  the  known  Baiker  codes  are  listed  in  Table  7. 


Table  7 

The  Known  Barker  Codes 

Length  of  code 

Code  Elements 

1 

+ 

2 

+  -  or  +  + 

3 

+  +  -  or  +  -  + 

4 

+  +  -  +  or  +  +  +  - 

5 

+  +  +  -  + 

7 

+  +  +  --  +  - 

11 

13 

+  +  +  +  +  --  +  +  -  +  -  + 
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From  Table  7,  it  is  seen  that  Baiker  codes  are  binary  codes  with  different 
lengths.  As  the  radar  is  made  operational,  one  code  length  is  selected.  Accor¬ 
dingly,  the  transmitted  wave  form  from  the  radar  is  made  of  segments,  the 
number  of  which  is  equal  to  that  of  the  bits  in  the  selected  code.  Here,  the 
“+”  and  signs  of  Barker  code  correspond  to  0-deg  and  180-deg  phase 
angles,  respectively.  They  are  used  to  set  the  phase  of  each  corresponding 
segment  from  the  transmitted  wave  form.  Take,  for  example,  the  application 
of  Barker  code  of  length  4,  i.e.,  +  +  -  +.  The  phases  of  the  first,  second  and 
the  fourth  segments  of  the  transmitted  signal  are  set  to  0  deg,  and  the  phase  of 
the  third  one  is  set  to  180  deg.  Thus  Barker  codes  are  phase  codes.  An 
onboard  microprocessor  is  used  to  select  the  appropriate  code  length. 

Let  the  duration  of  an  uncoded  transmitted  radar  signal  be  T.  Then  the 
range  resolution,  6^,  of  this  signal  is  (Eaves  and  Reedy  1987) 

s  =  c  ^  cT  (120) 

'■  2B  2 


where  c  is  the  speed  of  light  and  B  is  the  bandwidth  of  the  transmitted  signal. 

It  equals  to  1/T.  The  small  value  of  5^  implies  that  the  radar  can  resolve  small 
range  values;  therefore,  it  is  a  high-resolution  radar.  As  it  is  seen  from 
Equation  120,  a  high-resolution  radar  requires  that  B  be  very  large,  i.e.,  that  T 
be  very  small,  e.g.,  t.  This  requirement  leads  to  hardware  complications.  This 
important  criterion  can,  however,  be  achieved  through  pulse  compression  tech¬ 
niques  where  the  hardware  implementation  is  much  easier.  In  pulse  compres¬ 
sion  the  total  duration  of  the  uncoded  transmitted  radar  signal  is  divided  into  N 
pulses.  The  duration  of  each  pulse  is  selected  to  be  t.  Accordingly,  the  pulse 
compression  ratio,  CR,  is 


(121) 


and  the  radar’s  range  resolution  is 


(122) 


It  should  be  noted  that  the  values  of  the  bandwidth,  B,  in  Equations  120 
and  122  are  different:  it  is  larger  in  the  latter.  A  pulse  compression  technique 
is  implemented  in  the  design  of  the  high-resolution  radar.  In  this  case  the 
transmitted  signal  is  encoded  using  Barker  codes  as  described  previously  and 
thus  high  resolution  is  achieved.  Application  of  high-resolution  radar  includes 
object  detection,  target  classification,  train  mapping,  and  as  an  aid  in  any 
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application  in  distributed  clutter  suppression  (Eaves  and  Reedy  1987).  Two 
operations  are  performed  in  the  receiver  of  this  high-resolution  radar.  In  the 
first  operation  the  complex  reflection  coefficient  is  generated  as  a  ratio 
between  the  incident  and  reflected  energies.  In  the  second  operation  the  base¬ 
band  information  is  extracted  from  the  reflected  signal  through  envelope  detec¬ 
tion.  Information  from  the  first  operation  will  be  used  as  an  input  to  schemes 
used  to  solve  inverse  problems  under  study.  On  the  other  hand,  information 
generated  from  the  second  operation  will  show  the  changes  in  the  signal's 
amplitude  due  to  its  propagation  in  the  medium  under  testing. 


An  Improved  Portable  High-resolution  Radar 

In  this  section  an  improved  portable  high-resolution  radar  is  proposed. 
Some  information  regarding  this  radar  is  presented. 


Prior  art 

Up  to  the  present  time,  all  earners  used  in  radar  technology  have  been 
sinusoidal  signals.  However,  these  carriers  can  be  coded  or  not  coded  depend¬ 
ing  on  the  type  of  radar  used.  Codes  used  in  high-resolution  radar  are  Barker, 
Frank,  (^stas,  and  Welti,  among  others.  For  example,  carriers  of  ground- 
penetrating  radars  are  not  coded.  Most  ground-penetrating  radars  are 
somewhat  portable  but  not  to  the  extent  of  being  hand-held.  In  some  ground- 
penetrating  radars  the  antenna  is  designed  to  slide  on  the  surface  of  the 
medium  that  is  being  tested.  This  action  limits  their  applications.  As  another 
example,  a  synthetic  aperture  radar  is  a  high-resolution  radar,  however,  it  is 
not  a  hand-held  radar.  This  radar  takes  advantage  of  the  forward  motion  of 
the  airborne  radar  to  produce  the  equivalent  of  an  array  antenna  that  may  be 
thousands  of  f^t  long.  Moreover,  the  beamwidth  of  this  array  is  roughly  half 
that  of  a  real  array  of  the  same  length.  The  outputs  of  the  array  are  synthe¬ 
sized  in  a  signal  processor  from  the  returns  received  by  the  real  radar  antenna 
over  a  period  of  up  to  several  seconds  or  more. 


The  concept 

The  high-resolution  radar  system  is  designed  to  be  hand-held,  portable,  and 
microprocessor  controlled.  The  generated  carrier  is  an  exponentially  decaying 
sinusoidal  wave.  The  frequency  of  the  carrier  is  in  the  GHz  range.  In  this 
high-resolution  radar,  analog  to  digital  conversion  is  performed  through  an 
integrated  circuit  chip  at  a  rate  of  multiple  gigasamples  per  second.  In  addi¬ 
tion,  digital  signal  processes  are  applied  and  inversion  schemes  are  imbedded 
in  rnemory  chips  to  generate  electromagnetic  profiles  in  real  time.  In  order  to 
achieve  high  resolution,  digital  codes  such  as  Barker,  Welti,  or  Frank  codes 
cm  easily  be  implemented  in  this  radar  system  through  the  application  of 
C  MOS  logic.  The  coding  of  the  carrier  is  made  possible  through  the 
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application  of  digitally  controlled  phase  shifters.  Power  consumption  of  this 
high-rcsolution  radar  is  moderate.  Thus,  in  addition  to  being  hand-held,  the 
complete  system  is  mountable  on  any  vehicle  as  well.  As  a  re^t,  this  radar 
system  can  be  used  to  perform  target  and  media  identification  in  real  time. 

This  process  is  totally  performed  by  the  microprocessor  in  the  system,  and 
identification  results  are  displayed  visually. 

System  description 

The  portable  high-resolution  radar  is  of  four  parts:  a  transmitter,  a  receiver, 
a  microprocessor,  and  d-c  power  supply. 

a.  Transmitter.  A  block  diagram  of  the  transmitter  is  shown  in  Figure  1. 
The  transmitter  includes  the  following  items: 


(1)  A  keypad. 

(2)  Shift  registers. 

(3)  A  carrier  generating  circuit. 

(4)  A  phase  shifter. 

(5)  A  power  amplifier. 

(6)  A  transmitting  antenna. 


Figure!.  Transmitter 
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b.  Receiver.  A  block  diagram  of  the  receiver  is  shown  in  Figure  2.  The 
receiver  includes  the  following  items: 

(1)  A  receiving  antenna. 

(2)  An  envelope  detector. 

(3)  A  circuit  to  generate  the  complex  coefficient  of  reflection 
function. 

(4)  A  circuit  to  perfomi  inversion  and  generate  electromagnetic 
profiles. 

(5)  An  analog-to-digital  convertor. 

(6)  A  digital  signal-processing  circuit. 


Figure  2.  Receiver 


System  operation 

The  exponentially  decaying  carrier  is  generated  in  the  transmitter  and 
encoded  by  a  code  previously  selected  by  the  user.  The  code  is  input  through 
the  keypad,  which  is  a  part  of  the  front  panel.  As  the  carrier  is  generated,  com- 
m^ds  from  the  microprocessor  are  issued  to  the  digitally  controlled  phase 
shifter  to  encode  the  carrier.  Accordingly,  for  every  bit  of  the  code  the  car¬ 
rier’s  phase  is  changed.  The  duration  in  time  of  every  bit  of  the  code  is  much 
smaller  than  the  duration  of  the  transmitted  carrier.  The  encoded  carrier  is 
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amplified  by  a  power  amplifier  and  fed  to  a  transmitting  antenna.  A  reflected 
signal  is  receiv^  by  the  receiving  antenna  and  three  major  operations  are 
performed  on  it: 

a.  Envelope  detection  extracts  the  envelope  of  the  returned  carrier  and 
generates  a  baseband  signal  characteristic  of  the  target 

b.  The  returned  signal  is  utilized  to  generate  the  complex  reflection 
coefficient  function.  This  function  is  then  used  as  an  input  to  a  subcir¬ 
cuit  in  which  an  inversion  is  performed.  As  a  result  electromagnetic 
profiles  of  the  medium  such  as  permittivity,  conductivity,  permeability, 
and  susceptibility  are  generated. 

c.  The  returned  signal  is  digitized  and  digital  signal  processes  are  t^plied 
to  it  The  outputs  from  these  three  operations  are  features  of  the  med¬ 
ium  that  is  being  identified. 
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4  Conclusions  and 
Recommendations 


Conclusions 

Through  examples  presented  in  this  report,  it  is  demonstrated  that  the 
inverse  problem  is  an  ill-posed  problem.  In  two  formulations  of  this  problem 
discussed  in  this  report,  there  is  not  enough  information  avaUable  that  may 
lead  to  an  exact  inversion  procedure.  Thus,  approximate  inversion  procedures 
are  sought.  In  case  of  electromagnetic  propagation,  which  is  the  case  investi¬ 
gated  in  this  report,  solutions  of  inverse  problems  are  profiles  such  as  permitt¬ 
ivity,  conductivity,  and  permeability.  Such  solutions  have  applications  in 
target  recognition.  Specifically,  in  real-time  automatic  target  recognition,  the 
speed  by  which  these  profiles  are  generated  is  considered  important.  Further¬ 
more,  accuracies  of  these  results  have  vital  consequences  on  the  recognition 
process.  Deriving  a  solution  for  the  inverse  problem  in  dispersive  scattering 
theory  is  generally  difficult.  In  this  report  the  pseudospectral  Legendre  method 
was  used  to  find  the  dispersive  properties  of  a  homogeneous  medium,  i.e.,  the 
susceptibility  kernel  G(t),  by  means  of  scattering  methods.  The  method  is 
based  on  representing  the  solution,  i.e.,  G(t),  in  terms  of  its  value  at  Legendre- 
Gauss-Lobatto  nodes.  As  a  result,  the  inverse  problem  was  reduced  to  a  prob¬ 
lem  of  solving  a  system  of  algebraic  equations.  In  another  inverse  problem 
where  the  unknown  was  the  potential  function,  the  formulation  of  the  inverse 
problem  was  based  on  the  Gel'fand-Levitan  integral  equation.  However,  it  was 
decided  that  explicit  solutions  of  the  Gel'fand-Levitan  integral  equation  were 
seldom  easy  to  obtain.  Through  the  formulation  procedure,  the  solution  for  the 
potential  function  produces  several  valuable  and  useful  equalities.  Accord¬ 
ingly,  the  application  of  the  pseudospectral  Legendre  method  to  these  equali¬ 
ties  was  considered.  The  method  was  based  on  representing  a  kernel  l^(z,y)  in 
terms  of  its  value  at  Legendre-Gauss-Lobatto  nodes.  The  approximate  solution 
of  the  inverse  example  was  generated  under  conditions  produced  through  the 
differentiation  of  the  Gel'fand-Levitan  integral  equation.  As  a  result,  the 
inverse  problem  was  reduced  to  a  problem  of  solving  a  system  of  algebraic 
equations.  In  another  formulated  inverse  problem  the  solution  was  to  find  the 
scattering  potential  from  given  scattering  infonnation.  In  this  formulation  two 
Gel'fand-Levitan-Marchenko  integral  equations  were  associated  with  the 
Zakharov-Shabat  partial  differential  equations.  Two  kernels,  Aj(x,y)  and 
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i42(x,)'),  together  with  R(x),  which  is  the  Fourier  transform  of  the  complex 
reflection  coefficient,  were  used  to  express  these  integral  equations.  Again,  the 
pseudospectral  Legendre  method  was  used  to  solve  this  inverse  scattering 
problem.  The  method  uses  R(x)  as  an  input  and  generates  values  of  the 
unknown  potential  function  at  the  Legendre-Gauss-Lobatto  nodes.  Thus  the 
problem  was  reduced  to  one  of  solving  a  system  of  algebraic  equations.  The 
rapid  rate  of  convergence  of  the  pseudospectral  Legendre  approximations 
(GotUieb,  Hussaini,  and  Orszag  1984)  and  the  5  function  property  given  in 
Ch^ter  2,  “Application  of  the  Legendre-Gauss-Lobatto  collocation  method  to 
Zakharov-Shabat  equations,”  make  the  approach  used  to  solve  the  three  inverse 
problems  given  in  this  report  very  attractive.  The  method  is  efficient  and 
yields  accurate  results.  Numerical  examples  furnished  in  this  report  support 
these  claims. 


Recommendations 

The  new  Legendre-Gauss-Lobatto  inversion  scheme  formulated  in  this 
report  was  applied  to  solve  three  inverse  problems.  The  results  from  these 
three  solutions  were  compared  to  corresponding  solutions  from  other  methods. 
The  new  scheme  produced  accurate  results.  However,  the  unknowns  in  these 
problems  were  one-dimensional.  The  following  is  a  list  of  recommendations 
regarding  the  future  of  this  new  scheme: 

a.  Carry  out  a  study  of  the  effect  of  noise  on  the  performance  of  the  new 
scheme. 

b.  Perform  detailed  parametric  study. 

c.  Apply  the  new  scheme  to  the  two-dimensional  case  and  examine  its 
performance. 

d.  Develop  the  high-resolution  radar. 

e.  Utilize  the  high-resolution  radar  along  with  the  new  developed  scheme 
in  a  target  recognition  procedure  for  evaluation  purposes. 
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Appendix  A 
Notation 


B  Bandwidth  of  a  transmitted  signal 
c  Speed  of  light 
D(x,t)  Displacement  field 

Er(x,t),E^(x,t)  Left-  and  right-moving  electromagnetic  waves, 
respectively 

E(x,t)  Electromagnetic  field 

Ei/(x.t)  Second  partial  derivative  of  electrical  field  intensity 
with  respect  to  the  variable  of  t 

EjJx,t)  Second  partial  derivative  of  electrical  field  intensity 
with  respect  to  the  variable  of  x 

EJz,t)  Second  partial  derivative  of  electrical  field  intensity 
with  respect  to  the  variable  of  z 

G  Susceptibility  kernel 

j  Constant,  VT;  index 

L  Depth  of  medium 

Lj^t)  Legendre  polynomial  of  order  N 

L’f^t)  First  derivative  of  Lf^t) 

r(x,t)  Impulse  response  of  the  inhomogeneous  medium 
r(K)  Complex  reflection  coefficient 
R  Impulse  response 
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S  Scattering  matrix  of  transmission  line  junction 

t  Time;  variable  in  the  Legendre  polynomial 

T  Temporal  length  of  the  impulse  response;  duration  of  an 
uncoded  transmitted  radar  signal 

v(z)  Potential  fimction 

X  Depth  in  the  dielectric  inhomogeneous  medium 

and  Z{x+dx)  Wave  impedances  at  the  points  x  and  x  +  dx, 
respectively 

z  Depth  in  the  medium 

6^  Range  resolution  of  an  uncoded  transmitted  radar  signal 
Zq  Permittivity  of  free  space 
efzj  Dielectric  permittivity 
K  Wave  number 
Hq  Permeability 

prK;c;  Frequency  domain  reflection  coefficient  of  the  medium 

T  Duration  of  each  pulse  of  an  imcoded  transmitted  radar 
signal 

®  Radial  frequency 
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